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FOREWORD 


This  report  was  prepared  by  the  Wind  Tunnel  Branch,  Aircraft  Labora¬ 
tory,  Wright  Air  Development  Center,  under  Researoh  and.  Development  Order 
Dumber  !|5$~U29*  "Dimensional  and  Axially  Symmetric  Transonic  Flow." 

The  author  wishes  to  express  his  gratitude  to  Dr.  K.  G.  Guderley  for 
contributing  the  leading  theoretical  ideas  behind  this  study,  and  also  for 
his  continued  encouragement  and  guidance,. 


This  work  was  performed  while  the  author  was  assigned  as  project 
engineer.  Theoretical  Aerodynamics  Unit,  in  the  Wind  Tunnel  Branch,  WADC, 


ABSTRACT 


The  simulation  in  wind  tunnels  of  flows  with  a.  sonic  free  stream 
velocity  is  studied  from  the  theoretical  standpoint,  and  a  preliminary 
experimental  phase  is  also  described. 

The  solutions  which  are  valid  asymptotically  are  extended  to  in¬ 
clude  higher  order  terns  so  that  flow  patterns  at  distances  of  technical 
interest  from  bodies  writh  and  without  lift  can  be  discussed.  This  is 
done  with  a  view  toward  determining  the  essential  characteristics  which 
streamlines  near  a  wind  tunnel  wall  would  have  to  possess  if  the  simula¬ 
tion  of  Mach  One  flow  were  to  be  obtained,  and  also  toward  determining 
which  body  shape  parameters  are  needed  in  order  to  apply  the  correct  scale 
factors.  The  results  indicate  that  the  ratio  of  model  to  tunnel  sire 
v/here  a  single  parameter  is  sufficient  to  give  the  conditions  at  the  wall 
is  quite  large. 

A  set  of  preliminary  wind  tunnel  tests  was  made  with  Inclined  plane 
walls  in  the  test  section.  This  approximation  of  the  correct  streamline 
contour  yielded  data  wv ich  appears  to  be  in  good  agreement  with  the  theo¬ 
retical  results  at  '%ch  One. 

The  security  classification  of  the  title  of  this  report  is  UNCLAS3I- 


PU PLICATION  RSVIEY' 


This  report  has  been  reviewed  and  approved. 
FOR  THE  COMMANDING  GENERAL* 


R.  G.  TTCEGqwCr 
Colonel,  USAF 
Chief,  Aircraft  Laboratory 
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INTRODUCTION 


In  recent  years  the  interest  in  obtaining  data  for  models  operating 
in  th©  immediate  vicinity  of  the  sonic  velocity  has  been  quit©  presistent. 
Such  quantities  as  local  Mach  number  distribution  and  drag  have  been  shown 
to  have  a  smooth  behavior  as  l&ch  One  is  passed,,  but  it  is  not  clear  how 
the  values  of  the  lift  curve  slopes,  inter-action  coefficients,  duot 
efficiencies,  stability  derivatives,  and  other  characteristics  which  are 
sensitive  to  the  relative  sizes  of  the  subsonic  and  supersonic  regions 
behave  in  mixed  flows® 

Present  methods  utilising  free  flight  models,  transonic  bumps,  and 
very  small  models  in  large  wind  tunnels  have  provided  some  excellent 
qualitative  data,  but  there  does  not  seem  to  be  an  experimental  technique 
which  can  provide  the  degree  of  accuracy  normally  obtainable  from  wind 
tunnels  operating  outside  the  transonic  range*  However,  there  are  several 
techniques  currently  under  consideration  for  filling  this  gap. 

The  purpose  of  the  present  investigation  is  to  study  the  extent  to 
which  wind  tunnels  can  be  used  to  determine  the  aerodynamic  characteristics 
of  models  in  the  immediate  vicinity  of  Mach  One. 

The  approach  used  is  to  first  consider  the  distortion  of  the  tunnel 
walls  which  would  be  necessary  in  order  to  simulate  a  streamline  in  a 
free  flow  field.  In  the  type  of  theoretical  analysis  used  it  is  possible 
to  express  the  velocity  potential  as  a  •’basic"  potential  plus  a  series 
of  higher  order  particular  potentials  which  are  added  to  satisfy  the 
boundary  conditions  at  the  body®  The  basic  potential  is  th©  on©  which 
is  significant  at  large  distances  from  the  body,  and  therefor®  largely 
determines  the  wall  configuration* 

To  obtain  a  general  feeling  for  the  actual  effect  of  higher  order 
terms  due  to  the  shape  and  lift  of  bodies  on  th©  flow  fields  at  various 
distances  from  a  particular  model,  these  data,  have  bean  computed  for 
several  simple  bodies  and  the  results  are  compared  with  the  solutions 
given  by  the  basic  potential  alone® 


The  distances  from  the  bodies  where  the  influence  of  the  higher 
order  terms  become  of  the  same  general  proportions  as  the  experimental 
scatter  gives  a  fairly  good  indication  of  the  ratio  of  model  to  tunnel, 
size  where  tho  influence  of  the  tunnel,  wall  on  any  model  of  a  certain 
group  would  not  bo  significant  in  the  results e 

This  report,  which  represent e  the  first  port  of  a  larger  program, 
gives  an  account  of  the  analytical  study,  and  in  addition,  th®  results 
of  a  preliminary  experimental  program*  In  the  experimental  phase,  the 
ideal  wall  shape  was  approximated  by  a  nozzle  block  with  a  throat 
followed  by  a  long  straight  section  set  et  an  angle®  A  comparison  of 
these  results  with  the  analytical  data  seems  to  be  quite  encouraging 
despite  the  approximation  used. 
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SECTION  I 


ANALYTICAL  STUDY 


It  is  somewhat  more  convenient  to  treat  the  case  of  planar  flow 
first,  since  the  hodograph  differential  equations  are  linear*  The 
general  approaoh  is  that  used  by  Dr*  K*  G*  Guderley  (Reference  1)* 


Formulation  of  the  Problem 

The  differential  equation  for  the  velooity  potential,  ,  can  be 
simplified  in  the  oaae  of  small  perturbations  of  a  flow  with  sonio 
free  stream  velocity  toi 


(r+O^XX-^YY*0  (1) 

With  the  introduction  of  the  inclination  of  the  velocity  vector  Q  , 

3*  W 

and  the  dimensionless  perturbation  velocity  given  by  V'lr+i) ’{5.-U 
as  variables,  the  application  of  the  Legendre  transformation  yields  in 
the  hodograph  plane,  the  Tricomi  equations  0  (2) 

A  set  of  solutions  to  this  is  given  by* 


..here  n  =  m/2  and£  ■=  90^41?S  •  The  functions  fn  (  £  )  satisfy  the 

equations  f"C(H)+tl[($  "  g"''  fn=0 


(3) 


It  is  the  hypergeometric  function,  fj,  used  in  Reference  1*  The  values 
of  the  coefficients  are  to  be  determined. 

Working  with  these  functions  reveals  that  the  even  values  of  m  give 
warieuB  symmetrical  body  shapes.  The  addition  of  odd  val\:es  gives  asym¬ 
metrical  bodies* 


The  corresponding  points  in  the  physical  plane  are  given  by: 
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By  using  the  consistent  approximation  that  y  =  0  at  the  body,  it  is 
possible  with  these  expressions  to  prescribe  the  values  of  x  and  at 
M  points  on  a  given  body  from  a  resultant  system  of  (M  +  1)  equations 
in  (M  +  1)  unknowns  and  hence,  to  find  the  values  of  An  and  17  •  The' 

position  of  the  nose  of  the  body  is  given  by  the  branch  point  of  the 
streamlines  in  the  hodograph  plane.  Thi6  occurs  at  a  finite  value  of 
I1?!  ;  i.e.,  there  is  no  stagnation  point.  Therefore,  all  such  bodies 
have  cusped  leading  edges. 

It  should  be  mentioned  that  the  (*?  +  1)  expressions  are  non-linear 
and  futhermore  the  question  is  not  settled  as  to  whether  the  procedure 
converges  as  H  is  increased.  Consequently  it  is  possible  that  the  points 
at  the  boundary  where  matching  does  not  occur  will  oscillate  in  an  unde¬ 
sirable  manner. 

Streamline  Shapes 

The  vertical  deviation  of  a  streamline  from  its  original  position 

far  upstream  is  given  by  Ys  /  <£ydX  •  The  zero  subscript  denotes 

J-oa  5 

the  space  point  considered.  Because  of  the  differential  equation  (1), 
the  following  expression  can  be  written: 

ex  (r*  0  -^Yv]dxdY=o 

Using  Green’s  theorem  one  obtains 


dy  f  f  ydX“0 
C  2  yc 


If  a  closed  contour  is  chosen  which  is  given  by  the  lines 

Y.0,Y.V*.-00.  '"(s..  Sketch) 


and  the  line  AB  which  arises 
from  the  mapping  of  a  line 
£  =£0from  the  hodograph 
•to  the  physical  plane. 


All  terms  drop  out  except: 

CvLrC*#-- 


dY  +  f%]  dx 
dy  /.09  *•  jy*o 
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The  integral  on  the  left  gives  the  streamline  displacement  at  YsYq*£s£0 
If  we  transpose  to  the  hodograph  and  integrate 

*  % 

where  %(c0)  is  the  value  of  1)  for  Y*0  and  where  Yo  is  the  displacement 
of  the  streamline  iff  =  0.  In  symmetrical  flows, ^  gives  the  actual  tody 
shape*  It  can  readily  be  obtained  by  the  summation  of  from  the  nose 
to  the  value  of  P  considered. 

Gusped-  Body 

The  simplest  body  which  has  been  investigated  by  this  method  is 
given  by  the  potentials  for  n  =  -1  and  n  =  2.  This  is  a  cusped  profile 
with  a  constant  pressure  gradient  over  its  length.  For  such  a  body, 
the  transformed  velocity  potential  is  given  by 

<f>Hy+ 1)  ^a0w*[A-l  1?lf.+ A2i?2f2  ]  (6) 


and  the  stream  function  by:  (see  page£7  ,  Reference  1) 


-/3  'A  r  -S/l'  -p  1 

^  »-<!„ W*(y+l)  3 />*/£/  sion  (T)-8)lA.llylf‘^Ajrflfiji 

Along  the  body,  ^  -  0  and ■ equation  (7)  gives 


(7) 

(«) 


Substitution  into  (h)  yields  *  Aj>  /*•  Using  equation  (!}.),  (5),  and  '* 
(&)  the  pressure  fields,  deviational  ordinates,  and  body  profiles  were 
computed  for  Ap/A.j  =  10,  100  and  1,000.  These  patterns  were  compared 


with  the  patterns  given  by  "basic" 
A 2/^-1  =  100  are  shown  in  Figure  1 


potential  (n  =  -1),  and  the  results 
•  It  can  be  seen  from  the  comparison 


curves  of  this  figure  that  if  ±.01  is  chosen  arbitrarily  as  the  experi¬ 
mental  scatter  in  Cp  ,  if  (y  =  1  •!{.),  and  if  g00A_,  the  differences 

between  the  two  curves  will  be  of  the  order  of  the  scatter. 


*  If  Xe  is  the  length  of  a  body  from  the  nose  to  the  sonic  point  on 
the  body,  and  y$r,  the  body  ordinate  at  that  point,  and  if  the  value  of 
A^j  is  fixed,  the  bodies  -which  arise  from  various  values  of  Ag  all  have, 

the  same  value  of  Yo  fa  •  Accordingly,  there  exists  a  "bluntness 
factor";  i.e.,  blunt  bodies  which  have  a  smaller  actual  thickness  can  in¬ 


fluence  the  asymptotic  solution  to  the  same  extent  as  a  thicker  body  with 
a  higher  slenderness  ratio.  It  will  be  shown  later  that,  y  1/ki 
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is  a  aharaoteristio  parameter  whioh  will  be  the  same  for  all  affine  bodies 
which  hare  the  same  asymptotic  flow  fields* 

Wedge  Profile 

Quderley  and  Toshihara  have  computed  the  pressures  on  a  wedge 
at  Wash  One  for  both  the  symmetrical  (Referenoe  2)  and  asymmetrioal 
(Reference  5)  coses*  It  is  a  fairly  straightforward  procedure  to  find 
the  flew  field  in  the  vioinity  of  suoh  a  body  as  follows:  There  are 
three  eomrenient  forms  of  solutions  in  addition  to  the  set  of  solutions 
mentioned  in  the  previous  seotion*  They  are: 


hi" £JAJ 


9  (e-t»>ef/.v3j9) 


co 


hi*o 


9(1 


c  osfc*?*/) 


9 

6o 


(10) 


f) 


cos 


WTrJl 

9o 


(ii) 


where  f.t  is  the  same  hypergeometric  function  discussed  previously  and 
in  this  case  represents  singularities  spaced  along  the  line  If  m  0  at 


intervals  of  2  <f>  *  <fi t+<f)  j  satisfies  the  boundary  condition 

V'-Oat  Bm  ±Qm  ,  the  semi-angle  of  the  wedge.  In  <p^  and  the 

funotlon  is  given  by  3 ft. *>c > « (tfgl*) 

where  2i  is  a  linear  combination  of  Bessel  functions  of  order  1/% 

These  two  solutions fj)  ^and<£  are  used  to  satisfy  the  boundary  con¬ 
ditions  along  the  characteristic  in  the  supersonic  region  which  emanates 
from  the  point  •  0 tQm 

Again,  the  predominant  term  near  the  origin  of  the  plane  is 
£(m  »  o)  and  6ne  can  find  the  influence  of  the  other  terms  for  various 
values  of  and  Q  • 


This  too  was  carried  out  for  small  values  of  and  £  *  The  resultant 
flow  fields  are  shown  in  Figure  2a* 
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Comparison  of  th a  Two  Bodies 


It  is  of  interest  to  compare  these  two  bodies  whose  solutions  now 
have  been  derived* 

It  0O  is  the  ratio  Y eAo»  the  following  table  gives  the  numerical 
values  which  have  been  obtained  for  the  cusped  body  and  the  wedge t 


I4/3-I 
X5(6c)  a-1 

YS/X£  A-l 


Cusped  Wedge 

2.62  2.80 

17.95  22.00 


These  values  would  be  identical  if  the  ratio  of  the  values  of  the 
coefficients  is  given  by  ^-Iq/^-I^  “  1.070.  If  this  were  true,  for 

the  cusped  body  Xg  *  2.83  &nd  for  the  w®dg«»  X©  *  2.80  ©©"^^. 

These  values  are  quite  close  and  in  light  of  the  faot  that  the 
aocuraoy  of  the  numberioal  methods  used  was  about  ±1$,  it  can  be  said 
that  within  this  accuracy,  a  wedge  and  cusped  body  which  have  the  same 
lengths  and  heights  to  the  sonic  lines  have  the  same  effeot  on  the 
asymptotic  flow  field.  If  this  could  be  shown  to  be  the  general  case, 
the  problem  of  simulating  Mach  One  flows  would  be  greatly  simplified. 

Asymptotic  Solutions  in  the  Physical  Plane 


A  method  which  was  first  used  in  Reference  3  i®  quite  convenient 
for  considering  the  asymptotio  flow  fields  for  both  the  planar  symmetrio 
and  axisymmetric  cases.  The  general  technique  has  been  used  in  Reference 
It-  to  obtain  the  asymptotio  shock  wave  patterns  at  Mach  One. 

With  a  free  stream  of  sonic  velocity,  the  velocity  perturbation 
potential  differential  equation  in  the  physical  plane  can  be  written  ass 

(y+  1)  &XX  -  0yy  “  A  fly/3  •  0  (13) 

In  the  planar  case  A  ■  0  and  in  the  aadsymmetrio  oase  A  ■  1.  If  a  solo*, 
tion  is  assumed  of  the  form  ^  *  y3n“2  f(^  )  where  $  ■  (y'  ♦  1 jy**  the 
resultant  equation  is  f*  (n2£  2  -  f ’ )  _  n  |  f’ (5n-5+A)-f (3n-2)  (3n-3+A)  -  0 

It  has  been  determined  that  for  the  planar  oase  fi  m  and  for  the  aa&» 
symmetric  case?!  ■  I4/7  (f  is  found  by  numerical  integration).  Asymptotio- 
cally  the  shock  coincides  with  a  line  of  constant  J  •  In  the  planar  oase 
J  ■  2.018,  and  in  the  axisymmetric  oase,  /  «  2.2ljD.  The  streamline 
deviational  ordinates  are  in  the  axisymnetric  case  given  bys 

y-|y-5A[fZ-1|  (f  T  *|y-5/7  ?  (14) 

and  in  the  planar  oase 

y  -  5  y  J/5  jV"  |  /  < f  -2  /  i">j  -  5  yl/5?  (15) 
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Since  it  is  simpler  to  obtain  the  flow  field  at  a  given  value  of 
x  and  y  by  this  method,  it  was  used  extensively  in  this  study.  The  plots 
of  these  functions  are  shown  in  Figure  1. 

Similarity  Considerations 


If  two  separate  systems  are  related  by  the  following  transformation: 

.  t  --  -  -‘/a  - 

+  txyi  •  2,  Z,  $  (X/)  ,  X  •  Z,  X  ,Y  =  2,Z,  Y  (16) 


Equation  (13)  becomes 

9*  & 


?7Y“  A^y/y  =° 


The  points  in  the  two  systems  given  by  these  relations  are  called  "corres" 
ponding”  points. 

The  preservation  of  the  form  of  the  differential  equation  indicates 

that  if  <j|>  is  a  solution  for  a  flow  problem  in  the  X  ,  f  plane,  then  is 

a  solution  in  the  x,  y  plane.  IfZ^  f I, only  the  scales  of  the  two 

systems  are  different,  and  if  Z;s  f,  HjJ*  I  y  The  relations  given  above 
are  known  as  the  transonic  similarity  relations. 


Consider  the  case  where  two  different  affine  bodies  in  these  two 
flow  fields  are  producing  the  same  conditions  asymptotically;  i.e.,  if 
the  axes  of  the  two  systems  were  made  to  coincide,  the  velocities  at  a 
given  point  at  a  large  distance  would  also  coincide.  Then ^he  two  flows 
are  described  asymptotically  by  the  identical  forms  $  s  Y  *  * -f(t)  and 
—  _  ... 

£  T  Y  f  ( $1  .  This  may  happen  although  Z|,  and  Z^  are  not  equal 

to  one,  but  for  the  prescribed  similarity  to  exist  the  values  of  ^ and 

t  must  coincide  at  the  identical  point  in  space.  In  general  the 
"identical"  points  will  not  be  the  "corresponding"  points. 

Using  Equation  l6  and  the  asymptotic  expression  for  ,  the  follow- 

„  1/  “  (5^*2^  J  -t|  t 

ing  is  obtained:  =  Z  y  fe  Z  2)  Y  {[  _X  (  V  )  ) 

*  1  zz%  Z\v 


if  z 


z 


zJLf 


the  expression  reduces  to 


fffx'i'*  Jlnd 


opendent  of 


of  the  value  of  Zj  or  Zpr.  Hence,  with  this  relation  between  Zp  and  Zp, 

fba  values  of  and  and  aloe  v-  and  |»  coincide  at  the  identical  points Q 

5  *'  • 
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Now  oonsider  the  possibility  of  relating  the  dimensions  of  two  bodies 
which  give  asymtotically  the  same  flows  at  identical  points*  If  oorres- 


ponding  points  are  denoted  by  like  subscripts  3  x±  m  z  2  *i*  ^  2^*1  ^i 

and  if  ^  y  is  not  too  large  then  ydx  and  y^  c J<j>  -  dx®  Using 

yvcoKsT. 

li  n  »  3 

the  above  relation  between  T-.  and  Znr  7-?“™  ~  -7 

a  ^  Yi  "  y±Zi 


2  (1  «*  a) 


With  those  expressions  one  can  write 


().<  —  3)  + 

2(1  “  n) 


This  equation  ic  independent  of  Zy  if p  -  n/(ii  n  -  3)*  This  value  of 

/Ogives  the  power  law  relating  the  streamline  shapes  in  the  two  flows® 
The  relation  in  the  planar  case  (where  n  *  h/5)  is  given  by  ft  m  -ij.  and 
and.  in  the  axi symmetric  case  (where  n  *  h/l)  by  ft  I4/5 ® 


It  is  possible  to  go  further  and  to  say  that  if  the  values  of  x  and 
x  at  one  set  of  corresponding  points  are  denoted  by  the  subscript  1  and 

the  values  of  y  and  jr  are  taken  at  another  set  denoted  by  the  subscript 
.  a  _  —  $ 

2,  then  x±  yl>p  ”  yvj  •  Tiie  0686  where  point  1  is  at  the  nose  of  the 

profile  and  point  2  is  at  the  sonic  point  on  the  contour  corresponds  to 
that  discussed  previously  for  the  cusped  body. 

The  assumption  regarding  ^^y/^y  is  not  always  valid.  For  sufficient¬ 
ly  slender  bodies,  in  the  planar  case,  it  is  valid  everywhere,  and  it  it 
describes  the  way  in  which  bodies  of  an  affine  family  will  be  related  to 
give  the  same  asymptotic  flow  fields.  However,  in  the  axi  symmetric  case 
it  is  not  to  be  expected  that  this  assumption  will  hold  very  close  to 
singly  connected  body  shapes,  but  rather,  the  given  power  law  would  apply 
to  rings  or  stream  tubes  near  a  set  of  bodies* 


The  flow  field  in  the  vicinity  and  in  the  interior  of  closed  axi- 
symmetric  bodies  has  been  considered  previously*  Karman  (Reference  7) 

»  /.i  L 

a®  th®  pr@d £@siis8&B& 


gives  the  relation  ^  Jy 

expression  at  the  boundary*  The  results  ie^  «4|.  for  this  case© 


Correlation  of  B&nic  Planer  Solutions  in  the  Physical  Plane  with  Isdog: 


Solution® 

J.t  the  limitiur:  1C  oh 
hr  rtf?  solution  ~ "  1-  / 


e  tg  /?  'ff  V 
■  “■ 

r'~’~ 


*■  1  *  is 
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£•  (y*0 


••  •  Using  Equation  (i+)  for  x  and  y: 

^  X'fe  1  •(  •  A.,  (y+f)  f2  f£;+3f  £,)[-3(y*f1  %)^ 


For  9  •  Imh,  A-i  -  1.018.  Using  this  value  for  A-j.  gives  the  same  asymptotic 
nltMS  at  any  point  x,  y  in  either  system. 

Lot  Equation  (16)  be  applied  to  the  previously  discussed  cusped  body 

in  an  x,  y  plane  such  that  the  length  and  height  to  the  sonic  line  are 

given  by  Kg  “1.0  and  Y?  -  .1736.  This  body  corresponds  to  the  20° 
wedge  which  was  used  in  the  experimental  phase. 

The  flow  for  such  a  body  can  be  obtained  by  choosing  IgA-l  give 
the  proper  thickness  ratio  and  then  applying  a  scale  change  ,  Y*Z*Y 

2/3  1/6 

since  A^  ( A^/  )  (2.08)  and  Yg  =  A.^  (A^/  A ( )  (2.1+7) 

l/2 

15.  -1S_  .  *JI±6  .  M.  Z  47  or  jA?  /7.0 

xo  *3  1.0  A-I  *0*  A-l 

and  Z2  *  .0364 


If  a  streamline  of  interest  is  at  Y  =  3*00,  the  following  relations  can 


be  used  to  get  the  asymptotic  results  from  the  solutions  in  the  physical 

1/3  4/5 

plane,  X«Zj(y+l)  <Y/Z2)  •  1.6  5  £ 


,/5» 

Y-  5  Za(Y/Z*>  Y  • 


■$» 

44  1  Y 


-IAJ  . Z/5  •  .* 

Or*  *fy4l*  6f  /lt)  f  *  .2  5  6  f 

These  results  will  be  used  in  a  later  discussion. 

The  solution  given  for  the  basic  singularity  of  axisymmetric  flow 
should  at  sufficient  distances  give  the  correct  flow  fields  for  all  bodies 
of  finite  span.  Again  the  problem  is  to  determine  the  relation  between  the 
"strength"  of  the  basic  singularity  and  the  body  dimensions.  It  is  of 
interest  to  determine  whether  the  "sufficient"  distance  is  a  practical  one 
in  terms  of  model  and  tunnel  dimensions. 
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Lifting  Cases 


A  convenient  way  to  treat  bodies  at  small  angles  of  attack  is  to 
linearize  the  shift  of  the  mapping  in  the  hodograph  plane*  If  we  consider 
vhe  case  where  the  mapping  of  the  body  is  close  to  the  original  one,  and 

we  let  be  the  map  of  the  body  in  the  hodograph, 

lot  ^  *  Q  be  the  new  potential,  and  further  let  *  fjS)+F(Q) 

be  the  map  of  the  body  in  the  physical  plane,  Gtxderley  has  shown  in 

Reference  6  that  the  new  linearized  boundary  condition  is 

{I*  *  rja&0  *  F(0)-  F(fi) *  0 

if  the  body  location  or  contour  in  the  physical  plane  does  not  change* 

Cusped  Lifting  Bodies 

The  simplest  body  with  lift  that  can  be  treated  is  that  given  by  n  =  -1, 

-  •  at  •/, 

l/2,  and  2.  If  4>o  r  A-i 7  t,+V?  t,-(X6),  {1  «  A ,/tT)  fy  ...  (17) 

and  if  the  original  body  is  moved  through  an  angle  0<  ,  while  the  singularity 
remains  at  the  origin  of  the  plane,  the  boundary  condition  is  given 

t>y  {l-d^+^^g-O*  °  substituting  of  (16)  and  (17)  into  this 
expression  gives s 


It  is  possible  to  satisfy  this  equation  at  only  one  point  by  a  choice  of 
Al/2. 

/■  Reference  1,  page  54,  gives  recurrence  relations  between  4_  and 
H+  3/2  which  yield:  " 

6  ■  -A  V 

Using  these  relations,  if  Q  =  0  is  chosen  as  the  matching  point. 


Ayz  9  8«/& 
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V-(W0  <*•  A /A,) 


The  resultant  camber  shape  can  be  found  by  an  integration  of  the 
average  values  of  the  slopes  of  the  upper  and  lower  sides  of  the  body, 
and  the  numerical  results  as  shown  in  Figure  l\s.  indicate  it  to  be 
parabolic.  The  resultant  camber  line  shape  and  lift  data  are  shown 
in  Figure  !(.•  This  particular  case  may  not  be  too  enlightening,  but 
it  can  be  expected  that  the  chord  line  could  be  caused  to  pass  through 
any  desired  number  of  additional  points  on  a  straight  line  by  taking 
higher  order  terms  and  properly  choosing  the  coefficients.  The  resultant 
body  would  still  ha  Te  a  camber  line  which  is  not  straight,  but  the  de¬ 
viations  might  not  be  significant  in  the  results.  The  inadequacy  in  the 
nose  region  also  exists  in  the  other  approaches  to  the  lifting  problem 
in  transonic  flow  presently  in  use. 

In  regard  to  the  effect  on  the  walls  of  a  tunnel  on  lift  of  such 
bodies,  since  the  predominant  lifting  term  will  be  given  by  ^^and  Aft 

will  be  proportional  to  o(  ,  fortt>«  Go  it  is  apparent  that  the  lifting 
effect  will  be  of  lower  order  than  the  basic  singularity.  This  is  in 
contrast  to  the  subsonic  cases  where  the  influence  of  lift  predominates 
over  the  influence  of  thickness  at  far  distances. 

''"edge  Profile  with  Lift 


In  Reference  5*  the  method  of  linearizing  the  boundary  condition  is 
used  to  find  the  change  in  flow  due  to  angle  of  attack  on  the  previously 

determined  flow  of  a  double  wedge  airfoil.  In  this  case  the  ¥(Q)  -  0, 
i.e.,  the  wind  direction  is  changed  but  not  the  body  attitude.  The 
singular  point  then  moves  from  the  origin  of  plane  to  a  point 

y  =  0,  &  =  O'*  The  boundary  conditions  are  then  satisfied  along  the 

lines  0-  ±0.  and  also  along  the  characteristics  emanating  from  =  0, 

0- *0O-  It  is  shown  that  for  sufficiently  small  angles  of  attack, 

the  "separation  bubble"  which  must  occur  at  the  nose  has  negligible 
influence  on  the  lift. 

The  resultant  change  in  pressure  distribution  based  on  the  assumptions 
mentioned  is  shown  in  Figure  2c. 

With  these  results,  one  is  led  to  reason  that  a  "transonic  dip"  may 
exidt  in  the  dC^ycKycurve  for  wedges  with  thickness  ratios  greater  than 

about  .10.  Results  have  previously  been  obtained  analytically  for  the 
incompressible  lift  of  a  double  wedge.  Then  a  Prandtl-Glauert  correction 
is  applied,  this  gives  for  all  subsonic  Tach  numbers  higher  values  of 
dCj^/d  C(  than  those  obtained  at  I'ach  One.  (The  accuracy  of  these  results 

is  open  to  question  because  of  the  infinite  velocities  which  arise  at  the 
shoulder  and  at  the  nose).  The  values  obtained  for  purely  supersonic 
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flows  are  also  higher.  The  Mach  One  results  further  indicate  a  hori¬ 
zontal  tangent  to  this  curve.  Hence,  it  would  appear  as  though  the 
curve  has  a  rather  unusual  behavior  without  considerations  of  viscosity. 
This  is  one  of  the  problems  which  would  seem  to  justify  considerable 
effort  in  obtaining  experimental  data  at  Mach  numbers  very  close  to  one. 

In  most  practical  problems  with  three  dimensional  bodies  it  woulcl 
be  difficult  to  find  the  relation  between  the  basic  singularity  and  the 
body  dimensions  analytically.  In  light  of  this,  an  experimental  progra® 
has  been  initiated  to  help  fix  the  body  shape  parameters  which  would 
effect  the  flow  field  configurations  at  large  distances. 


SBC T I OH  II 


PRELIMINARY  EXPERIMENTAL  STUDY 


One  of  the  more  obvious  wayB  of  simulating  a  free  streamline  in  a 
wind  tunnel  is  to  have  flexible  '.'/alls  v/hic h  can  readily  be  adjusted  to 
a  desired  set  of  ordinates.  The  6“  supersonic  tunnel  at  ¥ADC  is  awaiting 
delivery  on  such  a  piece  of  equipment,  and  in  preparation  for  a  program 
utilizing  this  device,  some  preliminary  work  has  been  done  with  a  set  of 
fixed  geometry  blocks  which  have  a  throat  followed  by  a  long,  straight 
section.  See  Figure  5^>. 

The  initial  purpose  for  conducting  this  preliminary  study  was  to 
obtain  some  idea  of  the  influence  of  the  tunnel  boundary  layer  under 
the  pressure  rise  which  exists  ahead  of  a  model  in  Mach  One  flow.  Be¬ 
cause  of  this,  a  comparatively  large  model  was  used.  It  was  also 
desirable  to  determine  whether  or  not  it  was  possible  to  simulate  planar 
flow  by  shaping  the  portions  of  the  model  near  the  walls  so  that  the 
boundary  layer  on  the  •..all  would  not  change  the  effective  thickness  of 
the  model  in  that  region. 


It  will  be  noted  in  Figure  Ja  that  the  free  streamline  have  inflec¬ 
tion  points  quite  close  to  X  =  0.  Consequently,  one  might  obtain  a 
fairly  good  approximation  of  the  desired  flow  field  if  the  correct  slope 
is  simulated  by  the  straight  wall  at  the  model  abscissa.  It  is  necessary 
however,  to  evaluate  the  effect  of  small  changes  in  wall  contour  on  the 
flow  in  the  vicinity  of  a  model. 


I 


It  is  also  of  interest  to  determine  how  much  relatively  large  changes 
in  well  configuration  will  effect  the  flow  in  the  vicinity  of  a  given 
body.  It  can  be  seen  that  if  the  divergent  walls  which  were  used  in  this 
study  were  made  effectively  parallel,  the  result  would  be  a  "choked”  flow. 
This  corresponds  approximately  to  a  free  flow  at  some  Mach  number  which 
is  higher  than  the  upstream  Mach  number  in  the  choked  tunnel  but  less 
than  one.  It  is  further  evident  from  continuity  considerations  that 
with  parallel  -walls  and  a  symmetrical  pattern,  on  isobar  corresponding  to 
the  upstream  pressure  would  divide  thes e  other  isobars  v/r.ich  went  from 
the  wall  to  the  model  and  those  which  crossed  the  tunnel  centerline. 
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One  can  expect,  in  this  case,  that  if  the  pressure  drops  continuously 
on  a  model  toward  the  rear,  a  montonically  increasing  Mach  number  at 
least  to  the  sonic  line  along  the  wall  of  such  a  tunnel  will  exist. 

A  further  question  of  importance  involves  the  effect  of  lifting  • 
bodies  on  the  flow  patterns.  In  this  preliminary  phase,  the  wall  and 
model  pressure  distribution  were  obtained  for  a  20°  wedge  at  angles 
of  attack  up  to  6°. 

Experimental  Equipment 

1.  Wind  Tunnel 

This  tunnel  was  a  6"  square  test  section  and  is  powered  by  an  Allis 
Chalmers  1000  hp  axial  flow  compressor.  The  vertical  walls  of  the  test 
section  are  glass  panelled.  The  nozzle  blocks  are  equipped  with  static 
pressure  orifices.  The  block  ordinates  are  shown  in  Figure  5b.  Both 
mercury  and  TBE  manometer  panels  are  provided  for  obtaining  pressure  data. 

2.  Optical  Equipment 

The  system  for  flow  visualization  includes  interferometer,  schlieren, 
and  shadowgraph  with  arrangements  for  photographing  with  Polaroid  Land 
film  or  conventional  film.  A  comparison  of  interferograms  with  pressure 
data  indicates  that  at  the  density  level  used,  the  inherent  inaccuracy 
of  this  interferometer  is  about  ±.02  in  Faoh  number  and  consequently  it 
was  used  only  to  study  the  flow  pattern  qualitatively  in  this  program. 

Some  of  the  resultant  interferograms  and  schlieren  photographs  are  shovm 
in  Figures  8>,  9  and  12. 

3*  Wedge  Model  (See  Figure  10) 

A  10°  half  angle  single  wedge  of  1"  chord  and  5  1  /U"  span  was  used 
in  the  preliminary  runs.  In  an  attempt  to  simulate  planar  conditions 
as  well  as  possible,  the  3/8"  SaP  to  the  glass  windows  was  filled  with 
a  plastic  material  that  could  readily  be  shaped.  Fressure  taps  were 
located  at  mid-chord  at  positions  1/2"  from  each  tip  and  also  at  mid-span 
in  l/g"  intervals  chordwise.  The  assembly  was  mounted  to  a  rod,  which 
could  be  traversed  in  the  streamwise  direction.  The  shape  of  the  plastic 
fillers  which  gave  the  most  nearly  planar  flow  is  that  shovm  in  Figure  10. 
It  is  clear  that  perfect  simulation  of  planar  flow  cannot  be  obtained 
in  this  manner.  However,  the  difference  in  Cp  between  the  tip  orifices 

and  the  mid-span  orifices  never  exceeded  0.02  during  runs  where  data  is 

taken. 

No  after  body  was  placed  behind  the  wedge  section  since  it  was  de¬ 
sirable  to  minimize  any  influence  of  viscosity.  Angles  of  attack  were 
obtained  by  placing  shims  above  and  below  the  holding  screws  behind  the 
wedge  shoulder. 
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b»  Biconvex  Model  (See  Figure  12) 


A  biconvex  section  with  the  same  chord  and  thickness  as  the  wedge 
was  also  fabricated.  This  had  two  orifices  placed  at  the  mid-chord,  mid¬ 
span  position,  one  on  top  and  the  other  on  the  bottom.  The  same  plastic 
spacer  was  used  for  tests  with  this  model. 

Wind  Tunnel  Runs 


All  runs  were  made  with  PQ  =  1000  PSF  and  TQ  =  100° F.  Pressure 

data  were  taken  from  a  TBE  manometer.  The  following  is  a  listing  of 
the  configurations  of  the  runs  which  were  made  in  this  phase: 


Model 

in  Degrees 

Nose  Station 

Run  No. 

Wedge 

0 

I42"  to  5£" 

0  to  32 

wedge 

2,  6 

b2 "  to  5S" 

33  to  g2 

Wedge 

1/2,  1,  1  1/2 

b5 " 

S3  to  &9 

Biconvex 

0 

b2 "  to  5S" 

90  to  95 

Biconvex 

1/2,  1,  2 

I46" 

96  to  100 

Biconvex 

2 

b2n  to  5g" 

101  to  112 

Discussion  of  Experimental Results 


The  two  initial  problems  of  the  preliminary  phase  do  not  seem  to 
have  any  troublesome  consequences.  It  is  not  likely  that  separation  or 
appreciable  boundary  layer  thickening  will  occur  at  the  walls  of  a 
flexible  nozzle  in  the  regions  of  interest  even  with  comparatively  large 
models.  For  the  configurations  which  will  be  used  for  the  next  phase, 
the  simulation  of  planar  flow  by  shaping  the  plastic  fillers  at  the 
tips  is  quite  satisfactory. 

The  aerodynamic  data  obtained  appear  to  be  quite  interesting. 

Figure  6a  shows  the  wall  pressures  which  were  obtained  by  traversing 
the  wedge  model  from  Station  Lf2  to  Station  56  at  Of  =  O  while  Figure 
9a  shows  the  corresponding  schlieren  pictures,  ffith  the  ftodel  in  the 
most  downstream  position,  the  flow  is  choked  tjy  the  model  support  system. 
This  can  be  seen  from  the  fact  that  forward  movement  of  the  model  does 
not  influence  the  shock  position.  There  is  a  supersonic  region  down¬ 
stream  of  the  throat  which  is  terminated  by  a  normal  shock.  The  flow 
pattern  ahead  of  the  model  as  shown  in  the  schlieren  is  at  some  subsonic 
velocity  (M  =*  .75)  where  the  local  supersonic  region  first  begins  to  grow 
at  the  shoulder  of  the  wedge.  The  first  "choking"  by  the  model  is  evi¬ 
denced  by  the  forward  movement  of  the  tunnel  shock,  and  also  by  the 
appearance  of  sonic  pressure?  it  the  wall.  As  the  model  is  moved  forward 
the  supersonic  region  grows  until  it  fills  the  entire  rearward  field  of 
view. 
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It  will  bo  noted  that  the  pressure  distribution  at  the  wall  approaches 
the  theoretical  form  for  Mach  One  as  the  wedge  is  moved  forward.  At 
Station  i|6  for«  =  0  the  wedge  pressures  (Figure  5a)  as  well  as  the  wall 
pressures  (Figure  5°)  are  very  close  to  the  values  given  by  similarity 
theory. 

For  the  additional  tests  that  were  run,  this  position  was  considered 
to  be  the  optimum  configuration  which  cpuld  be  obtained  with  these  nozzle 
blocks.  It  can  be  seen  that  if  a  Mach  One  flow  can  be  sufficiently  well 
simulated  at  a  model  by  these  straight  blocks,  it  corresponds  to  a  free 
flow  over  that  model  with  some  other  body  or  "source  distribution"  far 
out  in  the  flow  field.  If  another  model  of  the  same  basic  singularity 
strength  is  substituted  for  the  original  one,  the  resultant  pressure 
distribution  at  the  wall  due  to  this  basic  singularity  should  not  change. 
Hence,  it  might  be  possible  to  find  the  body  shape  parameters  without  an 
ideal  wall  shape.  The  results  obtained  thus  far  have  been  quite  encourag¬ 
ing. 


In  Figure  5C  the  dotted  line  shows  the  pressure  distribution  given 
by  the  basic  singularity,  and  the  dashed  line  is  the  distribution  given 
by  a  cusped  body  of  the  same  length  and  thickness  to  the  sonic  line  as 
the  wedge  which  was  used.  Based  on  the  foregoing  discussion  of  the 
analytical  results,  the  theoretical  pressure  for  the  wedge  itself  should 
be  very  close  to  this. 

The  pressures  on  the  surface  of  the  wedge  at  Station  U6  are  in  very 
good  agreement  with  the  theoretical  values  of  Reference  3  for  at  least 
&0%  of  the  surface.  The  nose  region  is  not  accurately  described  by 
Reference  since  the  I*?  I  goes  to  infinity  under  the  approximations  used 
while  the  actual  flow  should  have  a  stagnation  point  at  the  nose. 

Dr.  Guderley  has  recently  corrected  these  results  to  include  the  nose 
region.*  The  corrected  pressure  distribution  is  also  shown  in  Figure  2b. 
The  agreement  with  the  uncorrected  theory  is  apparently  a  compensation 
of  the  approximation  errors  of  the  theory  and  the  thickening  effect  of 
the  boundary  layer. 

Wedge  at  Angle  of  Attack 

Figure  7b  shows  the  lift  data  obtained  by  recording  the  pressures 
on  the  surface  of  the  wedge  for  o(  from  0  to  6°.  Since  only  one  side 
of  the  wedge  had  enough  holes  to  obtain  the  pressure  distribution,  runs 
were  aetually  made  with  corresponding  positive  and  negative  angles  of 
attack.  To  check  whether  the  angles  were  the  same,  the  orifices  at  the 
mid-chord  point  of  the  opposite  side  were  also  recorded,  and  it  was 
necessary  that  the  pressure  at  this  point  match  the  curve  for  the  corres¬ 
ponding  run. 

It  will  be  noted  that  there  are  two  distinct  regions  of  this  lift 
curve.  The  dashed  line  is  the  theoretical  value  obtained  in  Reference 
5  for  £0  •  It  can  be  expected  that  as  the  supersonic  region  at  the 
nose  grows,  the  actual  curve  will  deviate  from  the  one  obtained  analytically 
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by  neglecting  this  region*  When  the  entire  upper  side  beeoaes  super¬ 
sonic  the  body  can  be  expected  to  act  essentially  like  a  flat  plate* 

The  flat  plate  follows  the  relation  and  henoe,  the  steeper 

slope*  For  thinner  wedges  the  range  of  angles  of  attaoh  where 
is  less  and  consequently  the  transition  will  occur  sooner* 

Biconvex  Section 

Figure  11  shows  the  wall  pressure  data  obtained  with  the  model  nose 
at  Station  46*  This  appears  to  have  the  closest  to  correot  pressure 
pattern.  Comparison  of  this  pattern  with  that  given  by  the  wedge  reveals 
some  interesting  information*  The  maximum  Cp  for  the  biconvex  section 
is  *300  while  for  the  wedge  it  is  From  the  similarity  considerations 

previously  discussed,  for  the  wedge  the  scale  factor  Z*  is  *03614.*  If  it 
is  assumed  that  the  ratio  of  the  values  of  Cp  would  be  the  same  asymptotic¬ 
ally  since  Cp  *  2  (f  +  l)”*^  f  ( — 2L— for  the  biconvex  section,  the 

value  of  Z2  would  be  .026I4.  The  interferogram  for  this  run  (Figure  12) 
indicates  that  the  sonic  line  is  at  X5  =  .i^O*  The  "corresponding* 
cusped  body  or  wedge  of  this  length  to  the  sonic  line  has  a  height  to 
that  Xs  ■  *113*  The  geometry  of  the  biconvex  section  giveB  Y5  ■  *113* 

Hence,  there  is  agreement  in  the  role  of  YsAs  for  all  three  seotions 
thus  far  considered* 

Conclusions 

A.  The  analytical  results  indicate 1 

1.  Within  the  accuracy  of  the  numerical  methods  used,  wedges  and 
cusped  bodies  of  the  same  length  from  the  nose  to  the  sonic  line  and 
thickness  at  the  sonic  line  give  the  same  asymptotic  flow  fields* 

2*  Planar  bodies  of  a  given  family  which  give  the  same  asymptotio 

4  . 

flow  fields  Imve  the  same  values  of  Yff  /^o  •  This  value  is  propor¬ 

tional  to  A.f  • 

3*  For  planar  flows,  over  bodies  of  normal  thickness  ratios^ the 
distance^ where  the  streamline  shape  will  be  given  sufficiently  well  by 
the  "basic"  solutionjis  within  the  bounds  of  experimental  feasibility* 

In  the  example  considered  it  is  about  20  ohord  lengths* 

ij.*  For  moderate  angles  of  attaok  the  influence  of  lift  on  wall 
configurations  is  small  compared  to  the  influence  of  thickness* 

3*  For  the  axisymmetric  case,  asymptotic  solutions  are  also  given, 
but  it  is  not  clear  how  to  correlate  these  results  with  particular  dosed 
body  shapes  except  by  experiment* 
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B.  The  Preliminary  Experimental  Study  Indicates: 

1#  Even  for  the  ratios  of  tunnel  height  to  model  thickness  less 
than  twenty,  which  were  used  in  this  study,  good  agreement  with  theo¬ 
retical  results  at  Mach  One  can  be  obtained. 

2.  The  straight  walled  tunnel  appears  worthy  of  further  study  as 
a  means  of  simulating  sonic  planar  flow  over  bodies. 

3.  With  due  consideration  for  the  crudeness  of  the  experimental 
setup,  pressure  distributions  on  a  wedge  witho^«  0o  are  in  good  agree¬ 
ment  with  the  predictions  of  Reference  3  and 

ij..  There  is  evidence  that  the  value  of  vs/xs  is  the  parameter 
which  will  determine  the  basic  source  strength  and  hence,  the  wall 
configuration. 


WADC  TR  52-gg 


16 


REFERENCES 


1.  Guderley,  K.  G.  Singularities  at  the  Sonic  'Velocity.  USAF  Tech¬ 
nical  Report  No.  F-TR-1171-UD°  19U2*  (Unclassified,  English) 

2.  Guderley,  K.  G.  and  Yoshihara,  H.  Flow  over  A  -"edge  Profile  at 
Mach  Number  One.  USAF  Technical  Report  No.  573^*  19^1$ •  (Un- 
classified,  English) 

3.  Guderley,  K.  G.  and  Yoshihara,  H.  Axial-Symmetric  Transonic  Flow 
Patterns .  USAF  Technical  Report  No.  5797.  September  19l+9»  • 
(Unclassified,  English) 

4.  Parish,  D.  T.  and  Guderley,  K.  G.  Asymptotic  Forms  of  Shock  haves 
in  Flows  over  Symmetrical  Bodies  at  latch  One.  •  USAF"  Technical  Report 
No.  6660,  March  1952*  (Unclassified,  English) 

5.  Guderley,  K.  G.  and  Yoshihara,  H.  Two  Dimensional  Unsymmetric  Flow 
Patterns  at  Mach  Number  One.  USAF  Technical  Report  No.  6653,  January 
1952.  ("unclassified,  English) 

6.  Guderley,  F.  G.  Two  Dimensional  Flow  Patterns  with  A  Free  Stream 

Mach  Number  Close  to  One.  USA1"  Technical  Report  No.  "  "T  ay  1951  • 

(Unclassified,  English)" 

7.  Karman,  Theodore  von.  The  Similarity  lav/  of  Transonic  Flow  Patterns. 
Journal  of  Mathematics  and  Physics.  October  19i|7»  "Volume  26, 

pp.  122-190.  (Unclassified,  English) 

2.  Perl,  and  Klein,  M.  Theoretical  Investigations  of  Transonic, 
Similarity  for  ^odics  of  Resolution.  NAG A  Technical  Note  No.  2?39, 
December  1950.  (Unclassified,  English) 


WADC  TR  52-22 


17 


FIGURE  1 
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b.  PRESSURE  DISTRIBUTION  (a«0,Y«l.4) 
e» «  IO# 


c.  CHANGE  IN  PRESSURE  DISTRIBUTION 
DUE  TO  ANGLE  OF  ATTACK  FOR  EACH  SIDE 
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FIGURE  2:  THEORETICAL  DATA  FOR  WEDGE  PROFILES 
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FIGURE  3 


D.  FUNCTIONS  FOR  AXISYMMETRIC  FLOW 


ASYMPTOTIC  SOLUTIONS  IN  THE  PHYSICAL  PLANE 
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FIGURE  i+:  THEORETICAL  DATA  FOR  CUSPED  BODY  T7ITH  LIFT 
A2/A-1  »  100  A  /k-i  -2  /  3 
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e.  WALL  PRESSURE  DISTRIBUTION 

<r 


FIGURE  5:  COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  DATA  FOR  UEDGE  FOR 
"WEDGE  AT  STATION  I46  (  -0) 
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a*o* 

FIGURE  6s  EXPERIMENTAL  TALL  PRESSURE  DISTRIBUTIONS  FOR  TRAVERSES  OF  20°  WEDGE 
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FIGURE  7*  EXPERIMENTAL  WALL  PRESSURES  FOR  20®  WEDGE  WITH  LIFT  (UPPER  WALL) 


FIGURE  9*  SCHLIEREN  PHOTOGRAPHS  OF  MODEL  TRAVERSES 
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